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CU MM ARY 

.Accurate calculation* to evaluate tho performance of the 
stoichiometric carb m-oxygen propellant system have been car- 
ried out for nozzle flow with and without chemical reactiono and 
with and without vibrational adjustment. The calculations show 
that, for frozen chemical flow, a lag of vibrational energy states 
at chamber conditions nearly doubles the reduction in lap, as 
compared v/ith flow in which complete vibrational equilibrium is 
maintained. On the other hand, lags in vibrational adjustment 
have practically no effect on the theoretical performance of hot 
propellant systems if chemical equilibrium is maintained during 
nozzle flow. The preceding conclusions are in agreement v/ith 
the results on other propellant systems obtained previously by 
use of an approximate evaluation procedure. 



-iii 



PART 


TABLE OF CONTENTS 
TITLE 


PAGE 




ACKNOWLEDGMENT 


i 




SUMMARY 


ii 


I. 


Introduction 


1 


II. 


Basic Equations 


3 


III. 


Outline of Calculations 


9 




A. Equilibrium Flow 


O 




B. Frozen Flow 


n 

/ 




C. Chemically- Equilibrium, Vibrationally 
Frozen Flow 


O 

/ 




D. Constant Composition and Vibrationally 
Frozen Flow 


14 


IV. 


Discussion of Results 


r. 




REFERENCES 


19 



- 1 - 



I. INTRODUCTION 

It is well-known that the overall performance of rocket pro- 
pellant systems is dependent on the extent to which physico-chemical 
changes occur during nozzle flow. Thu3 it is customary to report 
specific impulse data for "equilibrium flow" (the chemical reactions 
arc fast enough to maintain thermodynamic equilibrium at all times) 
and for "frozen flow" (the chemical composition remains unchanged 
at the concentrations corresponding to thermodynamic equilibrium 
in the combustion chamber). Practical experience has shown that 
this method of estimating reaconable upper and lower bounds for per- 
formance evaluation is generally in accord with experimental results. * 
Implicit in the customary propellant evaluation procedures is 
the assumption that complete thermodynamic equilibrium is main- 
tained during flow with respect to excitation of internal degrees of 
freedom. It i3 apparent, however, that sizeable departures from 
equilibrium internal excitation may occur, for example, for the vi- 
brational degrees of freedom, particularly for non-hydrogen contain- 
ing propellant systems in small motors for which the vibrational rc- 

(1 2 ) 

laxation times may not be short compared to the residence time. ' * ' 



*The usual simplifying assumptions in propellant evaluation, which 
do not involve consideration of rate processes, are the following: 
Thermodynamic equilibrium is reached in the combustion chamber 
after adiabatic reaction; expansion of the combustion products through 
the Laval nozzle i3 adiabatic; the products of combustion behave as 
ideal gases; the adiabatic expansion may be considered to involve 
one -dimensional flow (parallel to the nozzle axis) of non- viscous ideal 
gases; the velocity of the gases at the nozzle entrance position is neg- 
ligibly small compared to that at the exit position. For further details 
concerning customary evaluation techniques see, for example, "Quan- 
titative Evaluation of Rocket Propellants" by S. S. Penner, American 
Journal of Physics, Vol. 20, 1952, pp. 26-31. 
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\pproxi rate estimates of the effect of vibrational excitation on 

( 3 ) 

performance of rocxctc have been carried out previously. 

It is the purpose of the present analysis to present the first 
accurate evaluation of the probable combined effects of chemical 
and vibrational lags during nozzle flow. Performance calculations 
have been carried out for the following special cases: 

(a) Complete thermodynamic equilibrium is maintained 
("equilibrium flow"). 

(b) Complete thermodynamic equilibrium is maintained 
with respect to internal energy states but no chemical 
reactions occur during nozzle flow ("frozen flow”). 

(c) Complete thermodynamic equilibrium is maintained 
except that no adjustments of vibrational degrees of 
freedom occur during flow (chcmically-equilibrium 
and vibrationally -frozen flow). 

(d) Complete thermodynamic equilibrium is maintained 
except that no chemical reactions and no adjustment 
of the vibrational degrees of freedom occurs during 
flow (chemically- and vibrationally-frozen flow). 
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I he requisite basic equations for cases (a) to (d) have been 
(r ) 

given elcev/liere. * 1 or the present purposes it will be convenient 

to express the cherrical coupositi m by the rrole fractions for the 
various chemical species. In terms of mole fractions it ia readily 
shown* that equations (7), (3), (11) and (12) of reference (4) become, 
respectively: 

(a) For equilibrium flow 

n n 



1 — ? 

ir.' n“ 



' u“ = ill 
c e 



oc 

’oe 



(b) For frozen flow 



i IU“ = ah c 
Zee oc 



:=1 



n 






N? I I . -{51 /M )2Z nP H 

kc okc ' c e ksl kc o*vC 



(la) 



\ IJ . - il , 1 

kc L oac oke J 



(lb) 



(c) For cheinically-equilibrium and vibrationally-frozen flow 

i U u ? * = AH OC rs AH° C f (M /M ) 21. N? -HT lb l (lc) 

Zee e oe ' c e' k=l ke*“ ke kc - 1 ' ' 



(d) i or chemically- and vibrationally-frozen flow 



&A * * £ <, »«C 



oe 



n Uc’> 



( 1 < 1 ) 



In equations (la) to (Id), u represents the linear flnv velocity 
at the nozzle exit position, N® and are the equilibrium mole frac- 
tions of species k evaluated at the chamber exit and nozzle exit posi- 
tions, reflectively; 15 ^ and equal the total enthalpy (sensible 



or the relatively simple case of flow with frozen vibrational energy 
states, it is not necessary to use the analysis given in reference (4) 
since equations (la) to (id) can be written down immediately. 



plus chemical enthalpy) per oic oi pecie^. k if complete intcxnal 
equilibrium exists at the noszie entrance and exit positions, respec- 
tively; the superscript vib signifies that only the vibrational enthalpy 
is involved; the subscript o denotes complete energy equilibria 
exists for all constituents; the superscript o denotes cox. pletc ckea - 
ical equilibrium exists for all constituents. 

From equations (la) to (Id) it is apparent that the evaluation 

of u and hence of Isp * u /g requires the calculation of equilibrium 
c o 

mole fractions N? and N? as well as of total and vibrational en- 
ke kc 

thalpies. Since the exit temperature is determined by the assump- 
tion that the flow process is isentropic, it is obvious that the exit 
temperature T , and hence the numerical values of N? , , and 

C KC UC 

I^oke ’ depend upon the assumed type of flow process. 






in. Ui'i.m caloul r ions 

or the otoichion etric solid carbon-gaceous oxygen pro- 
pellant system mole of C{ c) reacting with 1 mole of the 

equilibrium mole fractions at any temperature may be determined 
conveniently by v/riting the overall chemical reaction in the lor n 

C + O (1-x) C0 2 4- x CO + [(x-y)/2] 0 ? 4- yO 

V/e arrive at equation (2) by assuming that the first chem- 
ical reaction taking place in the combustion chamber is 
which io followed by a dissociation of part of the carbon dioxide to 
carbon monoxide and oxygen according to the relation xCO. — *■ xC J 
+ (x/2) Og. This decomposition io then followed by a dissociation 
of part of the molecular oxygen to atomic oxygen according to the 
relation (y/2) O. - — *- yO. In equation (2) x is the number of moles 
of carbon diosdde dissociated and y/2 is the number of moles of 
molecular oxygen dissociated. 

From equation (2) we have 

n fc = (2 + x + y)/2 (3) 

where = total number of moles of reaction products. 

For the chemical reactions under discussion two equilibrium 
constants are obtained 

1 

K 1 = Po (p o 2 > 

‘2 = p CO < p o/ ^ P C0 2 ^ = K 9^ 10 



(5) 



art, in the notation of data co. npiled fro n 



where 'j, i\y and 

reference (5)* ami pj. is the partial pressure of species k. 

For an ideal gas 

P k = «i, P t /n, 

where is the number of moles of species k and p^ is the total 
pressure. 



( 6 ) 



By use of equation (6) in equation (4) we obtain 

1 



1 



Kj =y [(x-y)/2pT [2p t /(2 + x + yf (7) 

and, similarly, equation (5) becomes 



r 1 

K 2 = x|[(x-y)/2] 1 /(l-x)j [2p t /(2+x+y)] 



1 

7 . 



U) 



kquations (7) and (S) form a system of two equations in two 

unknowns, x and y. By solving the equations for the expression 

1 ' 

^P t (x-y)/(2 + x + y)~J ^ and equating the results it is readily shown 
that 



x 3 (p t -K2 + K x K 2 ) + x(3K 2 - KjE 2 )-2K 2 = 0, 



(9) 



and 



y = KjX 2 / 4- 21^ . (10) 

From equations (9) and (10), the quantities x and y, and the 
number of moles of each of the combustion products, can be deter- 
mined provided p. and T , the adiabatic flame temperature, are 

V C 

known. 



* Class notes for JP 200, California Institute of Technology, 1901. 

* * The effect of lag in vibrational energy states on the equilibria! 
constants has been neglected. It v/ill be shown in Section IIIC that 
the changes in equilibrium constant will not materially affect the 
numerical values obtained in the present calculations. 
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In the calculations for finding T the following relations are 



used 



Q av * " (n CO Ail £ CQ + n CO ? A I f co ^ n O AH f Q * 



O = 

C 



= V [ H okc ’ H ok(298. 16°K)] 



and 



Q uv = 



( 12 ) 

(13) 



where O is the heat released by the chemical reaction, Q is 
av 1 c 

the heat absorbed by the combustion products, and AH?^ ia the stand- 
ard heat of formation of species k. 



With equations (9) to (13) it is found by trial and error that 

the adiabatic flame temperature is 3750°It for p fc = 20 atmos. ^ 

The mole fractions as found by use of equations (9) and (10) 

and the average molecular weight M, as calculated from the relation 

M = ZN k M k (14) 

k. 

o 

are given at 3750 K and = 20 atmoe. and at various temperatures 
for p fc = 1 atmos. in Table I. Here Nj, and arc the mole fraction 
and molecular weight, respectively of species k. 

TABLE I. EQUILIBRIUM MOLE FRACTIONS N? AND AVERAGE 
MOLECULAR WEIGHTS M AS A FUNCTION OF PRES- 
SURE AND TEMPERATURE FOR THE STOICHIOMET- 
RIC CARBON-OXYGEN PROPELLANT SYSTEM. 



T,°K 


p fe , atmos. 


N c° ? 


ts CO 


% 


K o 


■vV 


3750 


20 


0. 2900 


C. 4442 


0. 1787 


0. 0u?f 


32. 3 




1 


0. 2502 


0.4631 


0. 1762 


o. no 


31. 39 


3000 


1 


0.4233 


0. 365 


0. 15 


0. 047 


’4. 91 


2750 


1 


0. 6190 


G. 2430 


1. 1150 


0. 01 6v 


30. 09 
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A. .quilibrium Flow 

The evaluation of the performance for equilibrium flow re- 
quires the introduction of an additional asoumption, namely, isen- 
tropic flow during the expansion procesD through the nozzle. This 
condition can be expressed in terms of the component entropies by 
the relation 



^<c (3 0 ke- RI "Pc- RW< tx )=(f ''c /: - ; C lf N Uc (S ok«- R1 “P c - RM, k e > < 15 > 

where 5^^ and S oke denote the total molar enthalpies at the cham- 
ber and exit conditions, respectively and R is the molar gas con- 
stant. 



The exit temperature, T R , can now be found by computing 

the right-hand side of equation (15) for several temperatures and 

comparing the results with the value of the left-hand side. By using 

the values of M and as given in Table I and the entropies taDu- 

lated in reference (5), the left-hand side of equation (15) is found to 

be 67. 65 cal. /°K-mole for T = 3750°K and p„ - 20 afcmos. and the 

c c 

right-hand side is found to be 66.43 cal. /°K-rnole and 68. 16 cal. /°A- 

mole at 2900°K and 3000°K respectively for p = 1 atmos. By linear 

interpolation it is found that T e = 2970°K. 

By use of equation (la) it is found that AH°^ = 22. COO kcal/ 

mole and 17.437 kcal/mole at 2900°K and 3000°K, respectively. 

oo 

Linear interpolation yields the result AlI Qe e 19. 04 kcal/mole at 

T s 2970°lv whence it follows that u s 2,200 meters /sec. From 
e e 




= u c /g 



the relation 



(16) 
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the specific impulse is seen to have a value -2 . 27 second* . 



1>. frozen Flow 

1 or frozen flow equation (15) simplifies to 
pSc< S ote- l ,7 » 

since = N° c and M g = M c> Qy proceeding as in Section UIA 
vc find 

Te = 2110°K; 

from equation (lb) 

AH C s 17. 12 kcal/mole 
oe 

is obtained; finally u g = 2105 meters /see. and Isp = 215 sec. 



C. Chemically-Equllibrium, Vibrationally-Frozen r low 

For chemically-equilibrium but vibrationaily-froaen flow, 
the isentropic relation which io used to determine T q becomes 



n - o » r 

ZI ( s 
k=l 



kc w okc - Klnp^Rl^) X (M c /M.) 

t K° [ S , -Rtap -RtaN° - (S^ b - ,:ri b )1 

_ ke L oke ke ' ke kc J 



(H) 



For diatomic molecules the vibrational contribution to the 

(7) 



entropy is, in good approximation, 
„vib 



= r|^ u £ exp(u) - lj - ln£l - exp(-u)J^ . 



(19) 



Here u = hcco*/kT, h represents Planck's constant, c io the velocity 
of light, k the Boltzmann constant, and <*> * the wave number corres- 
ponding to an energy transition from the ground vibrational energy 
level to the first excited vibrational state for the non-rotating diatom- 
ic molecule. For the linear triatoxnic molecule CO^, the total viora- 
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fcional entropy is £i en by the approximate relation 



(?) 



S vib = R 



i=l 



Uj ^exp(u.)-lj -In l-exp(-u i )]^| 



( 20 ) 



where u. corresponds to the four normal vibration frequencies of C0 2 
For the present calculations the numerical values of the cliaracterie- 
tic temperatures Q = hew* /k which are listed in Table II were u 3 ed. 



TABLE II. CHARACTERISTIC TEMPERATURES FOR 0 2 , CO, AND 



Molecule 


°2 ' 


CO 


C °2 


9, °K 


2228 


3067 


955.6, 955.6, 1988, 3493 



vib 

The thermodynamic function S , as found by use of equations (19) 
and (20), is listed in i able III for various temperatures. 



TABLE IU. S vlh (cal. /°K-mole) AS A FUNCTION OF TEMPERA 
TURE FOR O z , CO AND CC 2 / 





3750°K 


3000°K 


2900°K 


1600°K 


Com- 

ponent 


0,°K 


u 


r vib. 

S 


u 


^vib. 


u 


S vib - 


u 


^.vib. 

J 


°2 


2228 


. 5941 


3. 050 


0. 7427 


2. 628 


0. 7683 


2. 5597 


1. 393 


1.482 


CO 


3067 


.8178 


2.442 


1. 0223 


2. 027 


1. 0575 


1.9689 


1.917 


0.972 


1 C °2 


1988 


. 5301 


3. 271 


0. 6627 


2. 839 


0. 6055 


2. 7769 


l. 243 


1.677 


2 C °2 


955.6 


. 2548 


4.709 


0. 3185 


4. 270 


0. 3295 


4. 2025 


0. 597 


3.039 


3 C °2 


3493 


.9314 


2. 198 


1. 1645 


1.774 


1. 2045 


1. 7304 


2. 184 


0.769 


CC 2 Total 


--- 14.387 


--- 13.153 


--- 12.920 


--- 8. 5-14 
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Ay use of the apnropriate values of N^ and M listed in 
Table I and the appropriate values of S vl listed in Table III we 
can now proceed to find Te from equation (IB) by a trial and error 
method similar to the method described previously. With the val- 
ues already obtained we find the left-hand aide of equation (18) to 

be 67. 65 cal /°K-mole at T = 3750°K, and p =20 atmoo. The 

c c 

right-hand side of equation (18) is found to be ’67. 557 cal/°K-mole 
and 69. 012 cal/°K-mole at 2900°K and 3000°K, respectively. By 
linear interpolation the result T e = 29iO°K. is obtained. 

For diatomic molecules the vibrational contribution to 
the enthalpy is, in good approximation,^ ^ 

h -1 

H V1 = RTu £exp(u)-lj (21) 

For the linear triatomic molecule CO-, the total vibrational en- 
thalpy is given by the approximate relation: 

H Vib = KTZZ u. [ ejqptuj) -l] 1 (22) 

Using these relatione the values for H v ^ (in cal. /mole) listed in 
Table IV are obtained. 



TABLE IV. H ViJ AS A FUNCTION Of TEMPERATURE FOR 
0 2 , CO AND CO^. 





H v ^ (cal /mole) at 


Molecule 


3750°K 


3ooo°i: 


2900°K 


16CC°T 


°2 


5455. 8 


4C18. 0 


3829. 3 


1462. 3 




4815. 0 


3423. 6 


3242. 6 


104/. 1 


C °2 


23,247 


17,424 


16,869 


7116. G 
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has already been noted, t.ie effect of lag in the vibrational 

energy states on the equilibrium constants has been neglected. In 

order to see just ho v much of an cri’or this assumption introduces 

vve will now proceed with the computation of F^ for T = 3000°K by 

taking the effect of lag in the vibrational energy states into account. 

The starting relations for the reaction CO.,, — *CO + y O, at the temp 

eraturc T = 3O0O°K are 
e 

Al' e = - RTln X e (23) 

r 



and 



AF e = 



where v 



7* T lw. 16<°2>* al l98. 16< C0 > tiF «S. 16> C °2> * AT < re6ultin 8 

from the given chemical re 
action at 298. 16°K) (24) 

is the Gibbs free energy as defined by the relation 



F = H - TS, 



(2b) 



and the superscript e again denotes conditions at Coe exit temper- 
ature. The enthalpy, , and the entropy, S^ e , are found from 
the relations 



H. = H (H^ b - I-I^ h ) 
ke oke ' ke kc 



S, = S , - ( 

ke oke ' 



,vib „vib 



’ke 



kc 



) , 



(26) 

(27) 



respectively. By use of the enthalpies and entropies tabulated in 
reference (5) and of equations (25) to (27), together with the approp- 
riate values for S v * h and H vl ' listed in Tableo III and IV, the data 
listed in Table V are obtained. 



r 



l 3 y u : ys: i z - j,. 

, l II I MA P .? :rr ) : ' 

V = 3750'Y. 
c 





1 


2 






6 




ke 


r- 


kc 


“ F k(29S. 16) 


Y 


At ?.93. 16 


>n - 
nunent 


kcal/ 

mole 


cal/ 

V 


».ca.l / 
. >le 


kcal / 
ole 


-t.cal / 
mole 


kcal/ 
i. ole 


** 


44. 763 


; 1 . G' 


"’'•0. 22 


12. 99 


-94. 260 


+137. 23 


Co • 


25. 825 


65. 37 


171.74 


12. 03 


-32.308 


«• 1 


>•> 

Li 


26. 953 


60. 40 


17 \ 25 


12. 54 


0. 0 


-165.71 



In Table V, A’ ° represents. the standard free energy of formation, 
r ' n column (5) of 7.' able V we have 



AT (resulting fro. . the given chemical reaction at 29S. 16°K) 
= -3 ?. 603 - (-94. 260) = + 61. 432 



and, using equation (24), 



AF e = 4-6. 08 keal. fm ole. 



Using this value for AF e in equation (23) it is readily shown that 

A~ e 

hi iv- - - iU—B - 1. 018 

2 x 



or 

log K 2 = - 0.442 at ’ = 3000°!... 

3n the other hand, for complete energy equilibrium, 

log ,: 2 = -0.4692 at J? = 3000°K, 
log T- 2 a -0. 1091 at f = 3250°IC, 
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auc’ , ueivg an appropriate linear interp)l lim, 

log . ^ ^ - 0.44? at 1/ = . 3314 >: n * 3 or T = 3017°K. 



From the preceding result it it; evident t!iat neglecting the effect 

of lag in the vibrational energy states on the equilibrium constant 

K, at T = C0Gf° ' has introduced an error equivalent to 17°r. F.e>- 

ctition of the preceding calc* . 1 tionc sho./c that the error in 

will be s/cr. smaller. Therefore, we arrive at the conclusion 

that the error in the mole fractions at Te = 2910* ,» is small enough 

so that the resultant "e for tins erre i<? within 15° of the value whic 

would be obtained if t’^c effect f la^ in the vibrational energy states 

on the equilibrium constants bad been taken into account. 

Using the calculations given for case (a) it i9 found that 

AlI oc = +22. 26 kcal/mole at 1 - 2910°K from which it follows, by 

oe e 

using appropriate values of N° and El from Table I and of II V *‘ from 

Table IV in equation (lc) that i l A = AH OC = + 18. 72 kcal/molc, 

c e c 

u * 2,203 meters /6cc, and I = 225 seconda. 

6 C ^ 



D. Cons taut -Composition and Vibrationally-Frozen Flow 

“or constant -composition and vibrationally-frozen flow the 



isentropic relation, which is used to determine T , becomes 



n 



k=l 



n: 



kc 



[ 



S , - S , 

okc oke 



- < S kc - S ke > 



0 ■ 



Rln(pc /pe). 



(28) 



With the values of N^ c listed in Table I, the appropriate values of 
given in Table III, and the entropies tabulated in refei'cnce (j), 
it is found from equation (28), by use of a trial and error method. 
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thnt f c = 16 . 0 " 



sin.., the values if * j v buLited in '?able I, 
.vih 



tie appropriate value c of l v * " given in lable IV, and the enthalpies 

tabulated in refcn ice (b), in equatio i (Id) it is found th^.t i P u 2 = 

c» c c 

n 'i = 1,981 x/ictcro /*ec, and I n 202 seconds. 

e sp 
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iv. Drcub^ioN or results 

For convenience the results of the present analysis arc 
summarized in Table VI. Reference to Table VI shows that the 
effect of vibrationally-frozen flow is of no consequence when chem- 
ical equilibrium is maintained. In this case lag in vibrational ad- 
justment lowers the performance by roughly one percent. 



TABLt VI. SUMMARY Ol- R'-SULTo 



Type of T low 


T e (°X) 


u c (m/o) 


Isp^sec) 


Chemically-equilibrium and 
energy-equilibrium flow 


2970 


2220 


227 


Chemically-frozen and 
encrgy-equilibrimn flow 


2110 


2105 


215 


Che m ic a Uy - e quilibriu- r. and 
vibrationally-frozen flow 


2910 


2203 


225 


Chemically -frozen and 
vibrationally-frozen flow 


160C 


1981 


202 



On the other hand, the effect of vibrationally-frozen flow 
for frozen chemical flow is ouch as to decrease lop approximately 
6°/o. In this connection it is of interest to note that lags in vibra- 
tional adjustment, for a given propellant system, become more 
likely as the nozzle size is reduced because a reduction in nozzle 
size means increased cooling rates. In particular, it has been 
postulated that frozen-vibrational flow occurs in the small nozzles 
used in the pneumatic method for measuring temperatures in 
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combustion chambers . ^ 

The calculations on the carbon-oxygen propellant system 
have been repeated by using an approximate procedure, ^ which 
permits ready estimates of the effect of vibrational lag on per- 
formance for propellant systems for which the usual theoretical 
performance data are available. The results of these calculations 
are summarized in Table VII, together with data for other pro- 
pellant systems which were obtained by utilizing the concept of a 
constant average vibrational heat capacity during flow. ' 1 Refer- 
ence to Table VII shows that lag in vibrational energy states pro- 
duces a similar effect on performance for widely different propel- 
lant systems. A simple physical explanation for this observation 
is obtained by noting that frozen vibrational flow reduces the ef- 
fective heat capacity during expansion by nearly the same extent 
for chemically-frozen and for chemicaily-equilibrium flow. On 
the other hand, the effective heat capacity is much larger for flow 
in which chemical changes occur provided the propellant system 
is sufficiently hot. Hence the effect of vibrational lags will always 
be more pronounced for flow without chemical change than for 
chemicaily-equilibrium flow. Furthermore, the disparity will be 
largest for the hottest propellant systems. The preceding remarks 
are obviously in accord with the data given in T ible W which show 
a large effect of vibrational lag on performance for chemicaily- 
equilibrium flow only for the relatively cool RFNA-Aniline propel- 
lant system. 



